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Abstract 



In a type II superconductor the gap variation in the core of a vortex line induces a 
local charge modulation. Accounting for metallic screening, we determine the line 
charge of individual vortices and calculate the electric field distribution in the half 
space above a field penetrated superconductor. The resulting field is that of an 
atomic size dipole d ~ ea B z, a B = h 2 /me 2 is the Bohr radius, acting on a force 
microscope in the pico to femto Newton range. 
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The trapping of a magnetic flux $ Q = hc/2e by a vortex line in a type II superconductor is 
a well known phenomenon 0. Less familiar, however, is the fact that a vortex line in general 
traps an electric charge Q as well. It is the purpose of this letter to determine this vortex 
line charge quantitatively and to discuss the feasibility of its experimental observation. 

The vortex line charge Q has been discussed before by Khomskii and Freimuth and 
by Feigel'man et al. H (see also |Q) within the context of the sign change in the Hall 
coefficient, as observed in a number of type II superconductors |J. Here, we concentrate on 
the vortex charge and its accompanying electrostatic features, with a specific emphasis on 
its experimental observability. 

The main reason for the charge accumulation around the vortex is found in the particle- 
hole asymmetry, as quantified by the energy dependence of the density of states (DOS, per 
spin) N(E) at the Fermi level, Q cx dN(E) / dE\^. In the presence of particle-hole asymmetry, 
the carrier density n(/i, A) in the superconductor not only depends on the chemical potential 
fi, but on the energy gap A as well. The singular behavior of the phase at the center of 
a vortex leads to the formation of a core with a suppressed gap function A(R — > 0) — > 0, 
where R measures the distance from the phase singularity. With the chemical potential 
fixed, charge carriers (electrons/holes for a nearly filled/empty band) are expelled from 
this core region. Metallic screening drastically reduces the accumulated charge, however, 
in our analysis below, we show that the residual vortex line charge \Q\ ~ efc P (A TF /£) 2 is 
still experimentally observable (here, k F , A TF , and £ denote the Fermi wave number, the 
Thomas- Fermi screening length, and the coherence length, respectively; we define e > 0). 
In particular, one may envisage the classic geometry for the observation of vortices via 
the Bitter-decoration method ||, with the superconductor filling the half-space z < and 
penetrated by a magnetic field B || z, see Fig. 1. The vortex line-charge produces an electric 
field in the vacuum above the superconductor (z > 0), which corresponds to the one of a 
surface electric dipole d || z with unit ± charges separated by a distance ~ 1 A. The vortex 
charge is associated with the core size £ and therefore a much higher resolution can be 
expected in an electrostatic experiment as compared to the magnetic experiments probing 
structures on the scale of the penetration depth A L . In the following, we derive the vortex 
line charge Q and solve the "half-space" electrostatic problem for a single vortex and for the 
vortex lattice. Next, we discuss the observability of the vortex charge and close with a few 
more subtle questions regarding our analysis. 

The origin of the vortex charge can be understood on the basis of a textbook problem [|7| : 
Consider the Sommerfeld free electron model for a metal and determine the particle density 
n at fixed chemical potential /x, n(T) « n(0)[l + (it 2 /8)(T 2 / fi 2 )]. The density increase Sn is 
a consequence of the finite temperature smearing of the Fermi function combined with the 
finite slope in the density of states N' = dN(E)/dE\ fl (= 3n/8/i 2 in a 3D parabolic band). 
For a general Fermi surface with a smooth DOS we have 

Sn « (7rT) 2 Ar;/3. (1) 

Note that the sign of Sn depends on that of N' with Sn > for electron like carriers ||. 

Next, consider a BCS superconductor where the pair occupation probability v 2 . = (1 — 
£k/-Ek)/2 (£k — e k~ ^ an d E\l = (£k+^ 2 ) 1//2 denote the excitation energies in the normal and 
superconducting state, respectively) determines the density via n = 2J2k v l.- The opening 
of a gap A in the spectrum produces an analogous smearing in the occupation probability 
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and the density changes with A according to 



5n^ A 2 iV;in(^ D /T c ), (2) 

where u> D denotes the usual frequency cutoff on the attractive interaction. Indeed, the occu- 
pation probability in the superconductor resembles a Fermi distribution with a temperature 
T ~ A ||, in agreement with the results ([I]) and (0). 

In the presence of a vortex line the gap parameter turns to zero in the core, 

A 2 (R)^AlR 2 /(R 2 + e), (3) 

where R < A L denotes the radial distance from the phase singularity, Aoo is the magnitude 
of the gap parameter far away from the core, and the coherence length £ determines its 
spatial extent. The slow algebraic decay 5 A 2 ~ — A 2 ^ R) 2 is a consequence of the slow 
decay of the supercurrent j(R) ~ jo(,/R within the London screening length A L (j is the 
depairing current density; 5 A 2 drops to zero exponentially for R > A L ). 

We account for metallic screening within a Thomas-Fermi approximation, substituting 
/j, by the electrochemical potential ii + ep in the expression for the density n. The density 
modulation Sn(R) = n[ii + e<p(R), A(R)] — A^) is driven by the variation in the gap 
function A(R) and induces a scalar potential (fi(R), which is obtained from the solution of 
Poisson's equation 

VV(-R) = ^edn(R). (4) 
Linearizing in ip and A we arrive at 

[V 2 - \- 2 MR) = 4vre5n ext ( J R), (5) 
with the Thomas-Fermi length A TP = (87re 2 N fJj )~ 1 ^ 2 and the "external" density modulation 



ftw(fl) = -N A J^^ (6) 



(we have substituted the expression N'\n(huj /T c ), see @, by the phenomenological pa- 
rameter N^dliaTc/ dii using the BCS expression T c m hu> exp(— 1/N^V)). The integration 
of (§) over the planar coordinate R provides the total external line charge 

Q^^TreA^^ln^. (7) 

For a BCS model in the clean limit A^^N^ = kpfl/ii 4 with p, = m c ffV 2 /2 (for a nontrivial 
Fermi surface we have lm ^ jl in general), and using d In T c /d In fx « ln(hu) D /T c ) ~ 1 - 
10, we obtain an external line charge of order ek F , with only a weak dependence on the 
superconducting parameters T c and A L . 

We determine the real charge distribution p(R) = —e5n(R) = — V 2 <^(-R)/47r (positive for 
electron-like carriers) by solving the screened Poisson equation (||). In the limit A TP <C £ we 
obtain 

/ D x _ ea B d\nT c £ 2 - R 2 

P[R) ~ ^ih^w+W 3 ' { ) 
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where a B denotes the Bohr radius. Overall charge neutrality requires the total charge 
/ d 2 Rp(R) to vanish: The line charge accumulated within the vortex core is Q% = 
ea B (d In T c /d In //)/(47r£) 2 ~ Q ex tA 2 F /£ 2 and an equal and opposite charge is provided by 
the screening outside the core region, see Fig. 1. 

Next, we solve the electrostatic problem for a single vortex line penetrating a supercon- 
ductor filling the lower half space z < 0, see Fig. 1. The potential (p(R, z) generated by the 
charge density p ext (R, z) = —e6n ext (R)Q(—z) is obtained by solving the (screened) Poisson 
equation [V 2 — \~ 2 Q(— z)]<p(R, z) = —4iip ext (R, z). We decompose the potential into a bulk 
and an interface term, <p(R,z) = ip OQ (R)Q(—z) + (p (R,z), where (poo(R) denotes the bulk 
solution of (|5|). The interface term ifo(R, z) can be obtained from the Fourier Ansatz 

MR, z) = J 0^Vo (K) exp[*KR =F kfz), (9) 



with k+ = K and k~ = y K 2 + Atf referring to values z > above and z < below the 
vacuum - superconductor interface. Requiring continuity of <p and V<p> across the interface, 
we can express <£>o(K) through the source term tp^K) = ATrp ext (K)/(K 2 + A", 2 ). After 
transformation back to real space we arrive at the final expression 



p(R,z>0) = J —tp^K) 



d 2 K \K 2 + K 2 



xexp(iKR- Kz). (10) 

The integral is dominated by small wave vectors and we may neglect K 2 as compared to 
\~ 2 . Using p ext (K) w QextKo{K£)/ m (A L /^), with K the modified Bessel function, the 
integration over K yields 

m dliaT c min(z, A L ) 



e<p(R,z)[eV] w 0.8 ——4 In 



v [ J nil 

(R 2 + z 2 ) 3 / 2 ' { ' 

where we have chosen z > £ and all lengths are taken in Angstroms (note that (|ll|) is 
independent of the DOS, the latter appearing both in 5n ext and in the screening length 
A 2 P , and only weakly depends on the superconducting properties). The result ([11]) is the 
potential generated by a surface dipole d || z smeared on the scale £, (p(r) = d ■ r/r 3 , 

ea B z m dlnT c mxn(z, A L ) 
d=— -rr^ ln c • 12 

With the logarithms roughly compensating for the numerical 7r~ 2 , we find a dipole with unit 
± charges separated by ~ 1 A. The charge and field geometry are illustrated in Fig. 1. 

The corresponding electrostatic problem for a vortex lattice is solved in the same manner. 
The integral / d 2 K expfzKR — Kz] producing the dipole field has to be replaced by the sum 
over reciprocal lattice vectors K n of the vortex lattice, {2ir /a/\) 2 J2n exp[iK n R— K n z], where 
a A = (2/y / 3) 1,/2 ($o/-B) 1 ' /2 denotes the lattice constant. Usually the sum can be restricted 
to the six nearest neighbor lattice vectors and the result reads 
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m a 



e<p^(R,z)[eV] « 15.0 f-rr^ ln 

m e fi a A dmp 



1 d\nT c mm(z,a/\,\ 



X 



1 + exp(— 2irz/ciA) X! cos K n R 



(13) 



This completes our derivation of the charge- and the electrostatic field distribution for the 
individual vortex and the vortex lattice. 

Is this vortex charge observable in an experiment? The most straightforward attempt to 
identify the vortex charge is based on (scanning) force microscopy. Indeed, the observation of 
single charge carriers by force microscopy has been reported by Schonenberger and Alvarado 
ITJ . Below we consider two experimental setups: i) A grounded metallic tip is approached 



to an individual vortex. The vortex (surface) dipole induces a second dipole in the tip, 
resulting in a dipole - dipole attraction between the vortex and the tip. The expected force 
is estimated to be of the order of F in d ~ 1CT 17 N. ii) The metallic tip is biased against 
the superconductor. In this capacitor geometry, the bias voltage V drives a charge transfer 
from the superconductor to the tip, leading to a tip - surface attraction F ts which has to 
be compensated in the experiment. As the tip is approached to the vortex, the vortex- 
dipole - tip-charge interaction produces an additional force on the tip, which is the desired 
signal. The estimated force is proportional to the bias voltage V and is of the order of 
Fuas ~ 1CT 14 V[V] N. Note, that the mobile vortex dipole can be distinguished from static 
charged surface defects, as produced by adsorbed atoms and molecules, by driving the vortex 
with an external ac force and using a lock-in technique. 

In order to estimate the dipole-dipole attraction in the first setup i) we model the tip as 
a metallic sphere of radius p. Its center is chosen a distance ( > p right above the vortex, 
thus producing the maximal tip - vortex attraction. A straightforward calculation using the 



image charge technique (see, e. g., Ref. ||11|| ; we ignore higher order images) provides the 
result 

2n 2 + C 2 

F ind = 2pCrf 2 ^^- (14) 

Inserting the expression (|T2| ) for the vortex-dipole, using e 2 /a 2 i = 8.23 • 1CT 8 N, and choosing 
a typical geometry with p ~ (/2, we find 

F ind ~2.5-10- 8 (^ 4 N, (15) 

where we have assumed that {m/m c s){d In T c /d In p) ln[min(£, A L )/£] ~ 10. With p ~ 100 
A the resulting force is F m ^ ~ 2 ■ 10~ 17 N. 

Next, we consider the capacitor setup ii) where the tip is biased against the superconduc- 
tor. It is convenient to model this geometry with a tip of spheroidal shape. Using elliptic 
coordinates (see, e. g., Ref. |I2]]), we define the superconductor and tip surfaces through 
the coordinates r] sc = and r] tip = i] , respectively. Of the three parameters, the surface 
- tip distance (, the tip radius of curvature p, and the tip aperture 2$, only two can be 
freely chosen, p/( = tan 2 $ = (1 — i] 2 )/r] 2 . Solving Poisson's equation AV^ = 0, we find 
V{rf) = Vg(r])/g(r] ), resulting in an electric field E tip (a, R) = —2Vz/g(r] )(a 2 + R 2 ) 1 / 2 at 
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the superconductor - vacuum interface. Here, g(rj) = ln(l + rj)/{\ — rj) and a = (/r] is the 
scale factor in the transformation to elliptic coordinates. The energy of the vortex-dipole in 
the electric field of the biased tip is given by U(a,R) = — d ■ E tip (a, R)/2 (only half-space) 
and for the forces perpendicular and parallel to the surface we obtain 



(^bias'^bias) ~ ~7~T p2 /i I „2 / P2V3/2 ( 16 ) 




(see Ref. JL0J for the description of an ac technique used to separate the small modulation 
Fj^ ias due to the vortex-dipole from the large base force F ts due to the image charge attracting 
the tip). The above results apply for distances > £; upon approaching the tip closer to the 
vortex the details of the charge distribution become relevant and the forces change, e. g., in 
-^bias we have to replace the scale a by the spread of the vortex charge £. Using ( ~ p and 
$ = Ti /A we obtain the numerical expression 



Fbias~10- 9 V[V](^j N. (17) 

With p ~ 100 A the force amounts to -Fbi as ~ 3 • 10~ 14 V\V] N. At present, forces in the pN 
range are observed in state of the art AFM experiments and the fN regime is being attacked 
in the near future. 

A number of alternative experiments detecting the vortex line charge Q look promising 
as well. Here we only mention the basic ideas, i) One of the most sensitive electrometers 
is the single-electron-transistor (SET), e. g., see Ref. ||13|| , where the small central island 
is connected via tunnel junctions to the two leads. A capacitively coupled gate takes the 
device to its optimal working point. The device has to be fabricated onto the superconductor 
surface with only a thin insulating layer (~ 100 A) decoupling the two systems electronically. 
Vortices driven across the central island act to modulate the gate voltage via their line charge 
and the signal can be picked up via a lock-in technique. Using the setup i) above we can 
estimate the induced charge on the island to be ~ dp/( 2 ~ 10~ 2 e, which should be well 
detectable by an SET with a charge resolution of ~ 10 -4 e/-\/Hz |L3 . 
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Another straightforward idea is to imitate the original decoration technique of Trauble 
and Essmann ||, using electric rather than magnetic particles. It seems difficult, however, 
to imagine an electric analogue of the small ferromagnetic (Fe, Co, or Ni) particles, being 
assembled and spread onto the superconductor surface in a similar fashion. Alternatively, 
one may resort to the use of electrons, separated from the superconductor surface by a thin 
4 He film, and setting up a Wigner crystal. The interaction between the surface dipole array 
due to the vortex lattice with the 2D Wigner crystal will lead to new features affecting the 
physics of the vortex lattice and the Wigner crystal as well. 

We turn to some more subtle issues. On a phenomenological level we can express 
the external density modulation through the Ginzburg-Landau energy density T = 
iV M [a|A| 2 + §|A| 4 + 7|VA| 2 ], with the parameters a, j3, and 7 depending on the chemical 
potential /z, 5n CKt (R) = -S^l*- The term N^a)\A\ 2 oc ^ln(^ D /T c )(l - T/T c ) pro- 
duces the main contribution due to particle-hole asymmetry. A second term oc N' (1— T/T c ) 2 
originates from taking the derivative of the prefactor in T . Finally, a third contribu- 
tion from the gradient term, (A^//x)(l — T/T c ) |£VA| 2 oc (N^/p)(l — T/T c ) 2 , is present 



6 



even without particle-hole asymmetry and produces the well known London electric field 
E L oc WVg(R), compensating the centripetal force of the rotating superfluid. In our analysis 
above we have concentrated on the leading contribution in (1 — T/T c ) and in \n(hu D /T c ). 
Away from T c the other terms will slightly modify our result. 

In summary, we have determined the line charge associated with the formation of a vortex 
in a type II superconductor. The charge is mainly driven by the particle-hole asymmetry, 
its bare value is Q ext ~ ek F , and screening reduces this value down to Q ~ e/c F A^ p /£ 2 . We 
have solved the electrostatic problem related to the observation of the vortex charge on a 
superconductor surface and found the associated electric field to be that of a microscopic 
dipole d ~ ea B z. With our results we hope to motivate new experiments looking for the 
vortex charge itself (with an interesting relation to the sign change of the Hall effect), we 
propose a new imaging technique able to address the structure of vortex systems on the 
scale £, and we suggest to use the charge array set up by the vortices in other experimental 
areas such as the problem of electrons on 4 He films. 
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FIGURES 




FIG. 1. The superconductor (lower half-space) is penetrated by the magnetic field B. The 
resulting vortex line is charged due to the particle-hole asymmetry as quantified by the finite 
derivative of the density of states at the Fermi level, dN(E)/dE\ /Ji . Charge carriers are expelled 
from the vortex core (core radius ~ £, + region in the figure) and an equal and opposite screening 
charge on a scale min(aA,A L ) accounts for charge neutrality. The electric field generated in the 
upper half space (see field lines in the figure) is that of a surface dipole d of size ~ ea B . A tip 
approaching the surface is attracted to the dipole with a force depending on the specific setup, see 
Eqs. (15) for a grounded tip and (17) for a tip biased with a voltage V against the superconductor. 
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